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HYDRODYNAMIC  IMPACT  OF  A SYSTEM  WITH  A SINGLE  ELASTIC  MODE 
I— THEORY  AND  GENERALIZED  SOLUTION  WITH  AN  APPLICATION 

TO  AN  ELASTIC  AIRFRAME 1 
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SUMMARY 

Solutions  of  impact  of  a rigid  prismatic  Jloat  connected  by  a 
massless  spring  to  a rigid  upper  mass  are  presented.  The 
solutions  are  based  on  hydrodynamic  theory  which  has  been 
experimentally  confirmed  for  a rigid  structure. 

Equations  are  given  for  defining  the  spring  constant  and  the 
ratio  of  the  sprung  mass  to  the  lower  mass  so  that  the  two-mass 
system  proiides  representation  of  the  fundamental  mode  of  an 
airplane  wing.  The  forces  calculated  are  more  accurate  than 
the  forces  which  would  be  predicted  for  a rigid  airframe  since 
the  effect  of  the  fundamental  mode  on  the  hydrodynamic  force 
is  taken  into  account.  The  response  of  the  two-mass  system 
gives  the  response  of  the  represented  mode  and , although  no 
provision  is  made  for  taking  into  account  the  effect  of  secondary 
modes  on  the  hydrodynamic  force,  means  are  indicated  whereby 
the  results  may  be  used  to  approximate  the  response  of  modes 
other  than  the  fundamental  mode. 

Time  histories  of  the  hydrodynamic  force  and  structural 
response  are  given  for  wide  ranges  of  mass  distribution  and 
ratio  of  natural  period  to  the  period  of  the  impact.  By  use  of 
nondimensional  coefficients  these  results  are  made  applicable 
to  different  combinations  of  velocity,  weight,  angle  of  dead  rise, 
and  fluid  density.  Although  the  equations  permit  solutions  for 
different  combinations  of  flight-path  angle  and  trim,  an  approxi- 
mation is  given  for  correcting  the  results  for  the  combination  for 
which  solutions  are  given  to  other  conditions  within  a narrow 
range  indicated  to  be  of  primary  interest  to  the  design  engineer. 

In  a comparison  of  the  theoretical  data  with  data  for  a severe 
flight-test  landing  impact,  the  effect  of  the  fundamental  mode  on 
the  hydrodynamic  force  is  considered  and  response  data  are 
compared  with  experimental  data.  Consideration  of  the  funda- 
mental mode  alone  fails  to  account  for  the  fact  that  during  the 
impact  partial  failure  of  the  inboard-engine  mounts  occurred, 
but  use  of  the  theoretical  solutions  to  approximate  the  effects  of 
further  wing  torsion  leads  to  substantial  agreement. 

INTRODUCTION 

In  recent  years  the  development  of  large  airplanes  has 
caused  the  elastic  behavior  of  airframe  structures  during 
landing  impact  to  become  important.  The  work  which  has 
been  done  on  this  problem  has  been  handicapped  by  lack  of 
proper  knowledge  of  the  time  history  of  applied  ground 
reaction.  This  situation  has  been  particularly  acute  for 


seaplanes  because  of  difficulties  in  measuring  the  hydro- 
dynamic  force,  the  seaway,  and  the  manner  of  contact  with 
the  seaway. 

In  order  to  facilitate  the  interpretation  of  flight  data  and 
to  lead  to  the  prediction  of  design  loads  on  a rational  basis, 
a theoretical  hydrodynamic  study  was  made  and  tests  of  a 
rigid  float  were  conducted  at  the  Langley  impact  basin. 
Since  the  results  of  these  tests  agree  with  the  theoretical 
results  for  wide  ranges  of  the  pertinent  variables  in  numerous 
force  time  histories  (reference  1),  it  is  assumed  that  the 
theory  may  also  be  used  in  considering  the  effect  of  the 
upper-structure  elasticity  of  a seaplane  on  the  motion  and 
force  characteristics  of  the  hull  proper,  which  is  assumed  to 
be  rigid. 

The  bending  of  wings  during  impact,  which  for  modem 
flying  boats  is  the  primary  structural  action,  is  considered  in 
the  present  report  by  reducing  the  fundamental  mode  to  an 
equivalent  two-mass  system.  The  results  are  presented  in 
a form  suited  to  general  application  and  are  compared  with 
experimental  results  for  a particular  case.  The  equations 
showing  the  method  of  solution  are  included  in  appendix  A, 
and  a sample  data  sheet  is  given  as  table  I. 

SYMBOLS 

time  required  for  one-fourth  cycle  of  natural 
vibration 

t(  time  between  initial  contact  and  maximum  hydro- 
dynamic  force  for  rigid  structure 
t time  elapsed  after  initial  contact 
mL  lower,  or  hull,  mass  of  two-mass  system 
ms  upper,  or  sprung,  mass  of  two-mass  system 
m gross  mass  (T Vjg  or  msf-mf) 

W gross  weight 

g acceleration  due  to  gravity 

K spring  constant  of  spring  connecting  ms  and  mL, 

force  per  unit  deflection 

7i  i acceleration  normal  to  water  surface  of  nodal  point 
of  elastic  system,  multiples  of  acceleration  of 
gravity;  for  two-mass  system,  acceleration  of 
center  of  gravity 

na  oscillatory  acceleration  of  hull  about  center  of  gravity 
of  two-mass  system  or  nodal  point  of  represented 
mode,  multiples  of  the  acceleration  of  gravity 


1 Supersedes  N AC  A TN  1396,  "Solutions  Cor  Hydrodynamic  Impact  Force  and  Response  ol  a Two-Mass  System  with  an  Application  to  an  Elastic  Airframe"  by  Wilbur  L.  Mayo,  1847. 
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Vq  resultant  velocity  at  instant  of  contact. with  water 
surface 


P 

T 


7o 


0 


c, 

c, 

cd 

yanax 

/. 


mass  density  of  fluid 

angle,  of  trim;  angle  of  hull  keel  with  respect  to  plane 
of  water  surface 

flight-path  angle  at  contact;  angle  between  flight 
path  and  plane  of  water  surfaee 
angle  of  dead  rise 


nondimensional  time  coefficient  (t  V0  ) 

nondimensional  load-factor  coefficient  ^ j 


nondimensional  drafts  coefficient 

t 


(*J$D 


draft  at  instant  of  maximum  acceleration  . 
natural  bending  frequency 


Where  units  are  not  given,  any  consistent,  system  of  units 
may  be  used. 


THEORY 

HYDRODYNAMIC 

The  hydrodynamic  theory  used  in  the  present  report  is  the 
same  as  that  developed  in  references  1 and  2.  A basic 
differential  equation  which  gives  the  instantaneous  foree  in 
terms  of  the  instantaneous  position  and  motion  of  the  float 
is  given  in  reference  2.  This  equation  is  used  herein  to 
determine  the  effect  of  airframe  elasticity  in  altering  the 
motion  and  force  time  history  (appendix  A).  The  solution 
is  based  on  the  assumption  that  the  float  docs  not  change 
trim  during  impact.  In  this  connection  .the  pitching  moment 
may  be  large,  but  the  tune  of  the  impact  is  short  enough  to 
warrant  (at  the  present  stage)  negleet  of  the  resulting  angular 
velocities  and  displacements.  

The  solution  presented  herein  is  for  a prismatic  float 
with  sueli  beam  loading  that  the  chines  do  not  immerse, 
during  impact.  For  waves  that  give  the  severe  design 
condition  of  full-length  impact,  conventional  beam  loadings 
are  small  enough  to  cause  the  maximum  force  to  occur  at 
drafts  sufficiently  small  to  make  the  effects  of  finite  width 
and  chine  flare  secondary.  Reference  1 indicates  that  for  a 
conventional  float  neglect  of  the  pulled-up  bow  is  justified 
when  the  trim  is  3°  or  greater.  Although  for  high-trim 
landings  initial  contact  by  the  afterbody  may  substantially 
change  the  trim  before  the  . main  forebody  impact,  the 
negleet  of  afterbody  loads  is  justified  bocause,  during  the 
main  impact,  the  shielding  of  the  afterbody  by  the  forebody 
due  to  depth  of  the  step  and  to  keel  angle  is  such  as  to  mini- 
mize the  importance  of  afterbody  loads. 

STRUCTURAL 

A simplified  representation  of  primary  elasticity  of  an 
airframe  is  shown  in  figure  1.  A rigid  lower  mass  mL  is 
considered  to  be  connected  . by  a massless  spring  to  a rigid 
upper  mass  ms.  In  determining  the  fundamental  bending  of 
airplane  wings  part  of  the  wing  mass  must  be  included  in 
and  part  of  the  wing  lift  should  be  applied  to  %.  In 
the  present  report  the  gravity  foree  on.  each  mass  is  assumed 
to  be  balanced  by  wing  lift. 


Fiou&e  1 Simplified  representation  of  primary  elasticity  of  an  airframe'. 


The  problem  of  determining  the  properties  of  the  two- 
mass  system  so ' that  it.  is  representative  of  the  primary 
elastic  action  of  the  airplane  is  rather  simple  if  it  is  assumed 
that  during  the  impact  the  structure  defleets  with  the  shape 
of  its  fundamental  mode  of  vibration.  The  requirements  are: 

(1)  The  total  mass  of  the  simplified  system  must  equal 
the  total  mass  of  the  airplane  in  order  that  the  proper  nodal 
or  center-of-gravity  accelerations  can  he  obtained. 

(2)  The  energy  of  vibration  for  the  same  amplitude  of 
the  hull  and  lower  mass  (relativo  to  the  nodal  poinl)  must 
be  the  same  for  the  two-mass  system  as  for  the  considered 
mode  of  the  airplane  structure. 

(3)  Tho  natural  frequency  of  the  two-mass,  system  must 
be  the  same  as'  the  frequency  of  the  considered  mode  of 
airplane  vibration. 

Equations  which  permit  determination  of  the  masses  ami 
spring  constant  of  the  simplified  system  so  that  it  meets 
these  requirements  are  given  in  appendix  B.  These  equa- 
tions and  the  foregoing  requirements  aro  applicable  for  both 
landplanes  and  seaplanes. 

In  the  present  report  the  represented  structural  mode  is 
considered  to  bo  devoid  of  vibration  prior  to  the  instant  of 
impact.  Thus,  the  computations  may  represent  either  a 
first  impact  or  a subsequent  impact  resulting  from  a bounce 
sufficiently  high  to  cause  aerodynamic  and  structural 
damping  to  stop  the  vibration  during  the  time  the  seaplane 
is  in  the  air.  This  report  does  not  give  a representation  of 
successive  impacts,  such  as  might  occur  in  seaway,  which 
lead  to  accumulative  or  resonant  effects.  Available  flight 
data  indicate  that  a single  heavy  impact,  such  as  that 
considered  herein,  is  the  primary  cause  of  structural  failures. 

The  response  of  the  two-mass-  system  is  obtained  in  con- 
nection with  the  calculation  of  the  time  history  of  tho  hydro- 
dynamic force,  and  from  this  rosult  the  complete  resDonse 
of  the  represented  mode  can  be  obtained  by  the  simple 
procedure  given  in  appendix  B and  demonstrated  in  the 
section  entitled  “Comparison  with  Experiment.’’  The 
response  of  other  modes  to  the  force  computed  on  the  basis 
of  the  fundamental  mode  can  be  separately  determined 
and  superposed  (reference  3).  In  order  to  minimize  the 
complexity  of  the  solution,  however,  the  present  investi- 
gation does  not  provide  for  taking  into  aeeount  the  effect 
of  the  other  modes  on  the  hydrodynamic  force.  Although 
the  other  modes  may  have  a substantial  effect  on  tho  local 
loads  in  the  structure,  the  effect  of  these  modes  on  the 
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hydrodynamic  force  is  considered  to  be  secondary  as  com- 
pared with  the  effect  of  the  fundamental  mode. 

If  a large  number  of  solutions  for  the  two-mass  system 
have  been  made  in  order  to  determine  the  effect  of  the  funda- 
mental mode  of  different  wings  on  the  hydrodynamic  force, 
the  response  of  modes  other  than  the  fundamental  can  be 
approximated  from  the  use  of  a solution  for  a mass  ratio 
and  ratio  of  the  time  period  of  the  impact  force  to  the  natural 
period  of  the  two-mass  system  representative  of  the  con- 
sidered mode.  An  example  of  such  use  to  approximate  the 
effects  of  wing  torsion  is  given  in  the  present  investigation 
in  a comparison  of  results  of  computations  with  experi- 
mental results. 

If  the  response  of  more  than  one  mode  is  considered,  the 
structural  and  aerodynamic  damping,  which  are  not  con- 
sidered herein,  are  important  factors  in  determining  the  ex- 
tent to  which  the  maximum  response  of  the  different  modes 
should  be  superposed  without  regard  to  phase  relationship. 
It  is  expected  that  the  effect  of  the  damping  will  be  most  im- 
portant for  the  higher  modes  and  that  a result  leading  to  con- 
servative design  will  be  obtained  if  damping  is  not  considered 
and  the  maximums  of  the  first  two  or  three  modes  are  super- 
posed without  regard  to  phase  relationship. 

RESULTS 

Solutions  of  the  equations  in  appendix  A were  made  for 
wide  ranges  of  the  pertinent  variables.  Time  histories  of 
the  calculated  nodal  acceleration,  or  hydrodynamic  force  in 


terms  of  the  weight,  are  given  in  figures  2,  3,  4,  and  5 for 
ratios  of  the  sprung  mass  to  the  lower  mass  equal  to  0.25, 
0.60,  1.00,  and  1.36,  respectively.  Each  figure  is  three- 
dimensional;  the  third  dimension  is  tjtt  which  is  a ratio  of 
the  period  of  natural  vibration  to  the  speed  of  the  impact. 
In  representing  the  period  of  natural  vibration,  f»  is  taken  as 
the  time  required  for  one-fourth  of  a cycle.  The  speed  of 
the  impact  is  represented  by  making  tt  equal  to  the  time  be- 
tween initial  contact  and  maximum  acceleration  for  a'  rigid 
structure.  If  the  time  to  reach  maximum  force  for  the  elastic 
structure  should  be  used  in  defining  ff,  discontinuities  in  the 
time  to  reach  maximum  force  would  cause  discontinuities  in 
the  time-ratio  scales  of  the  plots.  (See  figs.  2 to  5.) 

An  expression  for  f,  may  be  obtained  from  the  relation 

and  equation  (A4)  in  appendix  A.  The  equation 


for  tt  is  as  follows: 


Ot,  /TTT\X/» 


(1) 


where 

time  coefficient  at  instant  of  maximum  acceleration 
for  rigid  body  (0.678  for  0=22.5°,  70=14°,  and 
t=3°) 

The  expressions  for  f,  and  tt  may  be  used  to  determine  that 


k 

U 


(2) 
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The  oscillatory  acceleration  is  given,  in  figures  6 to  9. 
These  figures  are  the  same  as  the  figures  giving  the  nodal 
acceleration  (figs.  2 to  5),  except  that  the  acceleration  plotted 
is  the  difference  between  the  hull  aceeleration.and  the  nodal 
acceleration..  The  time  history  of  the  hull  acceleration  can 
be  obtained  by  summing  the  two  plots  for  a particular  mass 
ratio.  On  the  basis  of  the  assumption  that  during  impact 
the  airframe  structure  deflects  in  a particular  mode,  a time 
Instory  of  the  acceleration  of  any  point  in  the  structure  is  j 
obtained  from  the  results  for  the  equivalent  two-mass  system 
by  the  following  procedure: 

(1)  From  the  deflection  curve  of  the  represented  mode, 
obtain  the  ratio  of  the  deflection  of  the  point  of  interest  to  the 
deflection  of  the  hull.  Both  deflections  are  taken  relative 
to  the  nodal  point. 

(2)  Multiply  this  ratio  by  the  oscillatoiy  acceleration 
given  cither  by  figure  6,  7,  8,  or  9 or  by  interpolation  between 
these  figures  for  the  mass  ratio  of  the  equivalent  two-mass 
system. 

(3)  Add  result  to  the  nodal  acceleration  given  by  figures 
2 to  5. 


Time  histories  of  the  acceleration  given  in  figures  2 to  9 
are  on  a nondimcnsional  basis.  The  nondinicnsionnl  co- 
efficients, which  contain  velocity,  weight,  fluid  density,  and 
acceleration  of  gravity,  were  used  in  reference  1 in.  a com- 
parison of  theoretical  data  with  impact  data  for  a float 
having  an  angle  of  dead  rise  of  22J°. 

Application  to  other  angles  of  dead  rise.— The  function  of 
the  angle  of  dead  rise  can  also  be  included  in.  the  nondiinen- 
sional  coefficients,  but  in  the  present  investigation  this 
function  is  isolated  and  treated  as  a factor  for  correcting  the 
results  presented  for  angle  of.  dead  rise  of  22J°  to  other  dead- 
rise  angles.  The  pertinent  relationships  between  results  for 
different  angles  of  dead  rise  may  be  expressed  as  follows: 


tcc  [mM)l 

(3) 

(4) 

where 

/(jS ) function  representing  variation  of  virtual  mass  for 

two-dimensional  flow  with  angle  of  dead  rise 
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f(A)  aspect-ratio  factor  including  effect  of  angle  of  dead 
rise  on  aspect  ratio;  ratio  of  virtual  mass  for 
three-dimensional  flow  to  virtual  mass  for  two- 
dimensional  flow 

np  acceleration  at  any  point,  either  oscillatory,  nodal, 
or  total 

The  shape  of  the  force  and  acceleration  curves  for  a given 
value  of  ijtt  is  independent  of  angle  of  dead  rise,  but  the 
effect  of  angle  of  dead  rise  on  as  given  by  relation  (3), 
does  enter  into  the  determination  of  the  value  of  tjtt  for  a 
particular  solution.  After  the  value  of  tjtt  for  a particular 
solution  has  been  determined,  the  acceleration  and  time 
values  for  an  angle  of  dead  rise  of  22^°  are  proportioned  by 
means  of  relations  (3)  and  (4)  to  the  corresponding  values 
for  the  angle  of  dead  rise  used  in  determining  f„ftt  in  order  to 
obtain  the  proper  acceleration  history.  . 

Although  adequate  impact  data  have  not  been  available 
for  checking  the  theoretical  equations  for  angles  of  dead  rise 
other  than  22}°,  the  theory  is  equally  applicable  to  planing 
floats.  Study  of  planing  data  has  shown  that  the  functions  of 


angles  of  dead  rise  used  in  equations  herein  are  approximately 
correct  for  angles  of  dead  rise  ranging  from  15°  to  30°^  Tfie" 
functions  are: 

Afl-Csy-1)’  ® 

Until  improved  functions  of  angles  of  dead  rise  are  obtained, 
functions  (5)  and  (6)  should  be  substituted  in  equations 
(1)  and  (2)  and  in  relation  (4)  to  correct  for  angles  of  dead 
rise  within  the  range  from  15°  to  30°.  Hough  approximation 
can  be  obtained  by  use  of  functions  (5)  and  (6)  for  angles  of 
dead  rise  greater  than  30°  but  not  for  angles  of  dead  rise  much 
less  than  15°.  In  reference  2 there  is  a discussion  of  the. 
inadequacy  of  f(A)  for  aspect  ratios  which  normally  occur 
for  small  dead-rise  angles.  For  very  large  angles  of  dead  rise 
and  moderate  velocity  the  static  forces,  which  are  not 
considered  in  the  present  investigation',  become  of  greater 
importance. 


272453-5' 
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Approximate  correction  to  other  flight  paths  and  trims. — 
The  flight-path  and  trim  conditions  of  primary  interest  to  the 
designer  tend  to  be  independent  of  the  variables  in  seaplane 
design  and  constitute  a limited  range  withiu  which  approx- 
imate correction  of  results  to  different. angles  of  flight  path 
and  trims  can  be  made  without  necessity  for  repeated  time- 
history  solutions.  For  any  particular  combination  of 
horizontal  speed,  rate  of  descent,  and  trim,  the  most  severe 
impact  load  for  most  of  the  structure  occurs  when  the  sea- 
way is  such  that  the  keel  contacts  a wave  slope  approximately 
parallel  to  it.  The  effective  angle  of  flight  path  and  trim  for 
such  an  impact  are  defined  relative  to  the  inclined  wave  slope; 
therefore,  the  trim  which  gives  maximum  force  is  zero.  The 
largest  flight-path  angle,  relative  to  the  keel  and  to  the 
critical  wave  slope,  is  also  associated  with  the  most  severe 
force.  The  value  of  the  largest  flight-path  angle  is  not  so 
definite  as  tlio  critical  trim. but  tends  to  bo  independent  of 
variations  in  size  and  wing  loading.  The  velocity  of  the  wave 
should  be  considered  in  determining  the  contact  speed  and 
flight-path  angle. 

The  equations  and  method  of  solution  given  in  appendix 
A permit  solution  for  different  flight-path  angles  and  trims; 
however,  approximate  correction  of.  the  results  in  figures  2 to 
9 to  other  positive  contact  angles  can  be  made  by  assuming 


that  the  proportionate  effect  of  the  structural  elasticity  on 
the  hydrodynamic  force  is  solely  dependent  on  the  ratio 
tn/tt.  Curves  given  in  reference  1 show  values  of  C,  f or 

aitun 

different  flight-path  angles  and  trims,  which  may  bo  sub- 
stituted in  equation  (2)  to  obtain  the  value  of  tjtt  for  differ- 
ent contact  angles.  In  making  the  approximate  correction, 
the  solution  presented  herein  for  the  obtained  value  of  the 
ratio  tjU  should  be  used  to  approximate  the  shape  of  the 
curve  giving  the  desired  time  history.  The  load  or  accelera- 
tion scale  should  be  corrected  to  the  different  contact  angles 
by  multiplying  the  present  result  by  the  ratio,  determined 
from  curves  given  in  both  references  1 and  2,  of  the  load- 
factor  coefficient  for  the  different  contact  condition  to  the 
load-factor  cocffieient  for  the  conditions  of  J3°  trim  and  14° 
flight-path  angle  considered  herein.  Correction  of  the  time 
scale  involves  a similar  procedure  in  which  the  time  coef- 
ficient is  used  rather  than  tho  load-factor  coefficient. 

Since  the  force  curves  for  a rigid  body  are  approximately 
the  same  shape  for  different  angles  of  flight  path  and  trims 
(reference  1),  the  approximate  method  of  correcting  to  dif- 
ferent angles  of  flight  path  and  trims  would  he  almost  correct 
if  the  structural  elasticity  did  not  affect  the  hydrodynamic- 
force  curve.  The . percentage  change  in  the  force  on  the 
float  due  to  elasticity  is  a function  of  the  percentage  change 
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in  draft  caused  by  the  elastic  compression  for  a given  eenter- 
of-gravity  position.  An  indication  of  tbe  validity  of  tbe 
approximate  correction  is  obtained  by  studying  tbe  extent 
to  which  tbe  ratio  of  tbe  spring  deflection  to  tbe  draft  is 
constant  for  impacts  of  tbe  same  values  of  tjtt  at  different 
angles  of  flight  path  and  trims.  From  the  expressions  for 
tjt{,  Ct,  Cit  and  C,  with  the  spring  deflection  assumed  to  be 
proportional  to  the  hydrodynamic  force,  this  ratio  may  be 
represented  by  tbe  expression  CiG,lGit  in  which  all  values 
are  for  tbe  instant  of  maximum  acceleration.  Tbe  variation 
of  this  expression  with  flight  path  for  an  impact  of  a rigid 
float  at  3°  and  12°  trim  is  given  in  figure  10.  Values  of  Ci 
and  C,  used  in  obtaining  this  figure  are  given  in  reference  1; 
values  of  Ct  were  obtained  in  conjunction  with  the  data  of 
reference  1 but  have  not  been  published. 

In  the  present  report  the  numerical  values  of  CiCtjCt 
have  no  significance  and  they  are  of  interest  only  because  of 
the  extent  to  which  they  are  constant.  Figure  10  indicates 
that  for  large  flight-path  angles  and  small  trims  the  ratio  is 
approximately  constant.  The  deviation  from  a constant 
value  of  this  ratio  is  due  to  planing  forces  which  exist  in  an 
oblique  impact  and  become  more  important  for  low  flight- 
path  angles  and  high  trims.  The  conditions  of  large  flight- 
path  angle  and  small  trim  previously  adjudged  to  be  of 
primary  interest  to  the  designer  constitute  the  ranges  in 
which  the  deflection  ratio  is  fairly  constant  and  the  approxi- 
mate method  of  correction  should  give  a fair  degree  of 
accuracy.  The  present  solutions  are  considered  to  be  for 
conditions  suited  to  correction  of  the  results  to  other  con- 
ditions of  greatest  practical  interest  ; they  represent  a moder- 
ately severe  combination  of  flight-path  angle,  wave  slope, 
and  trim  chosen  to  facilitate  correlation  of  the  theory  with 
an  impact  which  resulted  in  substantial  damage  to  a well- 
instrumented  flying  boat  during  flight  tests. 

The  equations  presented  herein  are  not  valid  for  zero 
trim;  an  assumption  that  the  float  is  prismatic  gives  solu- 
tions of  infinite  wetted  length  and  infinite  force  for  zero 
trim.  Solution  for  3°  trim  and  a prismatic  float  is  much 
simpler  than  a correct  solution  for  0°  trim  because  necessity 
for  consideration  of  bow  shape  is  eliminated.  Tbe  solution 
for  3°  trim  may  be  taken  as  an  approximation  of  the  critical 
design  load  or,  as  illustrated  in  tbe  following  section  in  a 
comparison  of  theory  with  experiment,  an  empirical  factor, 
which  includes  bow  effects,  may  be  used  to  convert  values 
of  acceleration  and  time  for  3°  trim  to  values  for  0°  trim. 

COMPARISON  WITH  EXPERIMENT 

Impact  basin  data  have  not  been  obtained  for  suitable 
models,  and  most  flight  landing  data  have  been  inadequate 
for  the  present  study.  The  only  data  which  appeared  suit- 
able for  this  comparison  are  those  which  were  obtained  with 
a four-engine  flying  boat,  the  data  for  which  have  not  been 
published.  Data  were  obtained  for  a large  number  of  test 
landings,  but  only  one  of  the  impacts  is  very  well  suited  to 
the  present  analysis.  This  impact  gave  loads  sufficiently 
high  to  cause  large  effects  of  elasticity  of  the  wings.  The 
impact  occurred  against  the  flank  of  a sizable  wave  (4  ft) 


Flight-path  angle,  jj,  deg 

Fictrar  10.— Variation  of  foree-penetratkjn  ratio  at  maximum  force  with  flight-path  angta 

and  trim. 

i 

and  thus  facilitated  tbe  use  of  results  based  on  a planar  water 
surface.  The  present  comparison  of  theory  with  experiment 
will  be  restricted  to  this  impact  since  other  impacts  involved 
more  complex  contact  with  seaway  and  gave  less  force. 

The  horizontal  speed,  rate  of  descent,  and  trim  were  re- 
corded. A large  number  of  pressure  instruments  distributed 
in  the  hull  permitted  determination  of  the  water  surface 
relative  to  the  hull.  Data  recorded  by  these  instruments 
indicated  that  the  wave  slope  in  contact  with  the  hull  was 
approximately  planar,  that  the  trim  relative  to  the  wave 
slope  was  0°,  that  the  resultant  velocity,  considering  the  speed 
of  the  wave,  was  85  feet  per  second,  and  that  the  flight-path 
angle  relative  to  the  wave  slope  was  14°. 

Structural  data  available  for  the  test  flying  boat  are  not 
adequate  for  the  present  analysis.  Since  this  flying  boat  has 
the  same  number  of  engines  and  approximately  the  same 
gross  weight  and  horsepower  as  a landplane  for  which  a 
large  amount  of  structural  data  is  available,  assumption  is 
made  that  the  flying-boat  wing  has  the  same  mode  shape  and 
mass  distribution  as  the  wing  of  this  landplane.  Use  in 
equation  (B6)  in  appendix  B of  data  for  the  landplane  given 
in  reference  3 leads  to  the  following  mass  ratio  of  the  two- 
mass  system  representing  the  fundamental  mode: 

—=0.25 

mL 

Based  on  study  of  the  accelerations  at  several  points  in  the 
test  flying  boat  during  periods  of  relatively  free  vibration  in 
which  the  fundamental  wing  bending  mode  appeared  to  be 
predominant,  a natural  frequency  of  3.6  cycles  per  second 
was  selected  for  use  in  the  present  example.  The  funda- 
mental mode  frequency  of  the  wing  of  the  landplane  is  3.4 
cycles  per  second.  (See  reference  3.) 

Since  the  mass  ratio  is  equal  to  0.25,  figures  2 and  6 are 
used  to  approximate  the  action  of  the  fundamental  mode. 
Further,  substitution  of  values  for  conditions  for  this  impact 
in  equation  (2)  results  in  use  of  the  specific  time  history 
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given  for  ^-—1.2,  Substitution  of  tlie  contact  conditions  in 
U 

the  load  and  time  coefficients  fixes  the  load  and  time  scales 
in  an  absolute  sense.  Correction  from  3°  trim  to  0°  trim  may 
be  made  by  assuming  that  the  shape  of  the  time  history  is 
approximately  the  same  for  both  conditions.  Empirical 
correction  of  the  curves  from  3°  trim  to  0°  trim  may  be  seen 
from  reference  2 to  require  a 10-percent  reduction  in  the 
acceleration.  An  analysis-  of  data  obtained  at  the  Langley 
impact,  basin  for  impact  at  0°  trim  indicates  that  correction 
of  the  time  scale  from  3°  trim  to  0°  trim  requires  a 10-percent 
reduction  in  tlie  time  values. 

Hydrodynamic  force. — The  nodal-point  acceleration  nt, 
obtained  by  the  procedure  discussed  herein,  represents  the 
hydrodynamic  force,  in  multiples  of  the  weight,  applied  to 
the  flying  boat.  Since  tlie  experimental  data  do  not  provide 
measurement  of  the  hydrodynamic  force  as  such,  direct  com- 
parison of  the  theoretical  force-curve  results  with  experi- 
mental results  is  not  permitted.  Instead,  a comparison  of 
the  theoretical  response  of  the  structure  with  the  experi- 
mental response  is  necessary,  and,  if  the  agreement  is  ade- 
quate, it  may  be  concluded  that  both  the  hydrodynamic  and 
the  structural  actions  are  adequately  represented. 

Before  a study  of  the  rosponse  of  the  structure  is  made,  the 
theoretical  effect  of  the  response  on  the  hydrodynamic  force 
should  be  observed.  This  observation  is  made  by  comparing 
the  force  curve  obtained  for  a mass  ratio  of  0.25  with  the 
force  curve  for  a rigid  structure.  Both  curves  are  included 
in  figure  11.  Tlie  curve  for  the  case  of  a rigid  structure 

( — =0  ) was  obtained  from  reference  1 for  3°  trim  and  cor- 

\mL  J 

reeled  to  tlie  conditions  of  the  present  example  as  previously 
indicated.  Comparison  of  the  curves  for  mass  ratios  of  0 
and  0.25  shows  that  in  the  present  example  tlie  theoretical 
effect  of  the  structural  elasticity  on  the  maximum  hydro- 
dynamic  force  is  to  reduce  it  15  percent.  For  the  hypo- 
thetical condition  of  a concentrated  wing  mass  located  at  a 
point  in  each  semispan  of  a massless  wing  structure,  the 
conditions  of  the  present  example  would  give  theoretical 
reduction  in  the  maximum  hydrodynamic  force  due  to 
structural  elasticity  of  44  percent,  This  result  is  indicated 
by  comparison  of  the  maximum  of  the  curve  in  figure  1 1 for  a 
mass  ratio  of  1.36  with  the  maximum  of  the  curve  for  mass 
ratio  of  0. 

The  curves  in  figure  11  show  only  reduction  of  the  hydro- 
dynamic  force  because  of  structural  plasticity.  . It  should 
not  be  concluded,  however,  that  the  effect  is  always  in  this 
direction.  Figures  2 to  6 show  that  in  some  cases  the  hydro- 
dynamic  force  is  increased;  the  maximum  increase  which  was 
calculated  was  of  the  order  of  12  percent. 

Hull  acceleration. — By  combining  the  results  for  the  two- 
mass  case  with  approximation  of  the  pitching  action  on  the 
basis  of  a rigid  structure,  the  following  equation  may  be 
obtained  for  tlie  bull  acceleration  at  different  longitudinal 
stations: 


0 .04  .03  .18  .16  .ao 

Time,  sec 

Fiansi  11 —Variation  of  hydrodynamic  force  with  time.  r-0°;  I'i-SS  feet  per 

second;  It'-lO.OOO  pounds;  f„— 8.0  cycles  per  second. 

where 

n*  hull  acceleration 

a distance  from  station  to  center  of  gravity 

k radius  of  gyration  (12  ft  in  present,  example) 

l distance  from  resultant  hydrodynamic  force  to 

center  of  gravity. 

The  forebody  length  of  the  flying  boat  is  31.75  feet.  Be- 
cause of  bow  effects,  a length  of  25  feet  is  assumed  to  .have  a 
rectangular  loading  for  this  zero-trim  impact  and  the  result- 
ant force  is  located  12.5  feet  forward  of  the  step  which  leads 
to  a value  of  l equal  to  8 feet.  Use  of  the  foregoing  procedure 
to  calculate  timo  histories  of  the  acceleration  for  two  stations 
in  the  hull  at  which  accelerometers  were  located  gives  the 
curves  in  figure  12.  The  maximum  accelerations  recorded  at 
these  stations  are  also  shown;  agreement  with  the  computed 
maximum  acceleration  is  good.  The  full  experimental  time 
history  is  not  included  because  tlie  film  speed  was  not  great 
enough  to  permit  accurate  determination  of  the  shape  of  the 
time  history.  This  factor,  together  with  some  uncertainty 
in  defining  the  exact  instant  of  contact,  prevents  exact  check, 
of  the  time  to  reach  maximum  acceleration ; thus,  the  experi- 
mental points  in  figure*  12  arc  located  at  the  theoretical 
time  of  maximum  acceleration. 

Although  the  agreement  of  calculation  with  experiment  in 
figure  12  is  good,  this  agreement  can  be  interpreted  as  con- 
firmation of  the  elastic  action  of  the  structure  only  to  the 
extent  that  disagreement  of  hull  accelerations  computed  on 
the  basis  of  a rigid  structure  can  be  shown.  Values  of 
maximum  accelerations  computed  for  a rigid  structure  are 
included  in  figure  12.  Tlie  disagreement  with  experiment 
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Fegubk  12.— Variation  of  bow  and  step  acceleration  with  time. 


is  greater  than  for  computations  in  which  the  elasticity  of  the 
structure  is  considered.  The  difference  is  small,  however, 
because  the  effect  of  elasticity  in  increasing  the  hull  accelera- 
tion for  a given  force  is  largely  offset  by  the  effect  of  elasticity 
in  reducing  the  hydrodynamic  force.  Since  the  probable 
accuracy  of  available  data  must  be  considered,  experimental 
proof  of  the  theoretical  effect  of  the  elasticity  on  the  hull 
acceleration  cannot  be  claimed.  Support  for  a conclusion  in 
this  respect,  however,  can  be  obtained  by  comparing  the 
theoretical  response  of  the  wing  with  the  experimental 
response.  If  the  wing  responds  as  assumed,  the  basic  equa- 
tions require  that  the  hydrodynamic  force  and  hull  accelera- 
tions be  as  calculated. 

Elastic  axis. — Use  of  the  nodal  and  oscillatory  accelera- 
tions of  the  representative  two-mass  system  to  predict 
accelerations  along  the  elastic  axis  of  the  wing  requires 
consideration  of  the  fact  that  the  elastic  axis  of  the  wing  is 
not  at  the  center  of  gravity  of  the  flying  boat.  An  approxi- 
mate correction  may  be  obtained  by  multiplying  the  results 
for  the  two-mass  case  by  the  factor/  which  is  used  in  equation 
(7).  This  correction  is  not  entirely  consistent  with  that 
given  by  equation  (7),  but  each  approximation  is  considered 
more  accurate  for  its  particular  case.  An  improvement  to 
the  present  correction  which  would  change  the  results  2 to  3 
percent  might  be  made,  but  the  complication  is  not  considered 
to  be  warranted. 

In  the  present  example  the  value  of  / is  0.86.  Application 
of  this  factor,  of  equation  (B7)  of  appendix  B,  and  of  per- 
tinent structural  data  given  in  reference  3 gives  the  curves 
in  figure  13  in  representation  of  the  acceleration  time  histories 
of  the  elastic  axis  of  the  wing  for  the  hull,  inboard-engine, 
nodal,  outboard-engine,  and  tip  stations. 

Wing  torsion. — The  torsion  of  the  wing  during  impact  may 
have  substantial  effect  on  the  acceleration  of  engine  and 
nacelle  masses  forward  of  the  wing.  Use  of  the  procedure 
which  gave  the  acceleration  time  histories  in  figure  13  to 
calculate  the  acceleration  at  the  engine  gives  a maximum 


acceleration  of  3ff  at  the  inboard  engines  and  a maximum 
acceleration  of  5.6 g at  the  outboard  engines.  The  relative 
magnitude  of  these  values  is  in  strong  disagreement  with 
the  fact  that  during  this  impact  partial  failure  of  the  inboard- 
engine  mounts  occurred,  but  the  outboard-engine  mounts 
were  not  damaged. 

Since  the  structural  data  used  in  the  preceding  computa- 
tion are  for  the  actual  fundamental  mode,  their  use  in  impact 
calculations  involves  assumption  that  the  coupling  between 
the  torsion  and  bending  corresponds  to  the  coupling  which 
exists  in  natural  vibration.  Actually,  the  torsional  deflec- 
tions in  impact  are  determined  not  merely  by  the  bending 
deflections  but  also  by  the  large  nodal  acceleration,  which 
does  not  exist  in  natural  vibration. 

In  the  present  example  the  results  which  have  been  calcu- 
lated for  impact  of  a two-mass  elastic  system  will  be  used  to 
predict  the  response  of  the  engines.  The  procedure  for 
doing  this  is  to  select  the  proper  solution  and  then  to  adjust 
the  acceleration  and  time  scales  of  the  two-mass  solution  so, 
that  the  maximum  acceleration  of  the  nodal  point  corre- 
sponds to  the  maximum  acceleration  of  the  elastic  axis  at  the 
engine  station.  The  acceleration  of  the  upper  mass  of  the 
two-mass  system  then  represents  the  response  of  the  engines; 
however,  because  of  the  eccentricity  of  the  impact,  an  in- 
crement must  be  added. 

Data  obtained  from  Guggenheim  Aeronautical  Labora- 
tory, California  Institute  of  Technology  include  the  torsional 
deflection  of  a station  inboard  of  the  outboard  engine  for  a 
given  moment  applied  at  the  wing  tip.  For  a flying  boat 
the  absence  of  cut-out  for  the  landing  gear  tends  to  give  a 
stiffer  wing;  therefore,  in  the.  present  example,  deflection 
measured  on  the  wing  of  the  landplane  slightly  inboard  of 
the  outboard  engine  is  considered  applicable  to  the  outboard- 
engine  station.  Relative  deflections  between  the  inboard 
and  outboard  engines  are  estimated  as  follows: 

(1)  Torsional  deflection  at  inner  engine  equal  to  1 unit 
due  to  each  engine,  or  2 units  total 

(2)  Torsional  deflection  between  inboard  and  outboard 
engines  equal  to  2 units  because  of  greater  distance  of  flexure, 
boosted  to  3 units  because  of  increased  flexibility  of  structure. 

On  the  basis  of  data  from  California  Institute  of  Tech- 
nology the  average  of  the  static  moments  of  the  inboard  and 
outboard  engines  is  taken  as  22,000  foot-pounds.  After 
determination  of  the  torsional  deflection  at  the  outboard 
engine  for  this  average  moment  applied  at  the  wing  tip, 
multiplication  by  the  ratio  1/2  gives  a value  for  the  static 
deflection  of  the  inboard  engine  and  multiplication  by  the 
ratio  5/4  gives  a value  for  the  static  deflection  of  the  out- 
board engine. 

An  approximation  to  the  calculation  of  the  response  of  the 
engines  to  the  total  acceleration  of  the  elastic  axis  is  to 
neglect  dynamic  interaction  between  the  engines  and  treat 
each  engine  as  a single-mass  oscillator  having  a natural  fre- 
quency determined  by  its  static  deflection.  Such  a proce- 
dure gives  values  of  7.6  cycles  per  second  for  the  inboard 
engine  and  5.1  cycles  per  second  for  the  outboard  engine. 
Use  of  these  frequencies  to  determine  tn  and  division  of  f*  by 
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Figcbe  is.— Variation  of  accelerations  at  stations  along  tbe  wing  elastio  axis  with  time. 


ti  ( f<  equal  in  this  case  to  the  time  to  reach  maximum  acceler- 
ation of  the  curvos  in  figure  13  for  the  station  in  question) 
leads  to  valuos  of  tJU  for  the  inboard  and  outboard  engines. 
The  next  step  is  to  select  a mass  ratio  which  for  these  values 
of  tjtt  has  a shape  of  the  nodal  acceleration  time  history, 
which  approximates  the  shape  of  the  acceleration  time  history 
for  the  elastic  axis  at  the  station  in  question.  In  the  proserit 
example  the  mass  ratio  of  0.25  is  used.  Scale  factors  for  both 
the  load  and  time  scales  axe  determined  so  that  the  maximum 
nodal  acceleration  for  the  two-mass  solution  will  agree  with 
the  maximum  acceleration  and  time  to  reach  maximum 
acceleration  of  the  elastic  axis  at  the  station  in  question. 
After  these  factors  are  applied  to  both  the  nodal  and  oscilla- 
tory curves  for  the  selected  mass  ratio  and  time-period  ratio, 
use  of  the  results  and  equations  (B7),  (B5),  and  (B2)  in 
appendix  B to  calculate  the  acceleration  of  the  sprung  mass 
of  the  two-mass  system  gives  accelerations  of  the  engines. 
Approximation  and  superposition  of  the  pitching  action  on 


the  basis  of  a rigid  structure  leads  to  the  solid-line  curves 
given  in  figure  14  for  the  accelerations  of  the  engines. 

Also  included  in  figure  14  is  the  design  ultimate  acceleration 
for  the  engine  mounts.  Comparison  of  the  calculated  engine 
accelerations  with  this  value  shows  agreement  j>f  the  calcu- 
lation with  the  fact  that  partial  failure  of  the  inboard-engine 
mounts  occurred  but  the  outboard-engine  mounts  were  not 
damaged. 

During  the  impact  an  accelerometer  was  located  at  the 
outboard-engine  station  intermediate  between  the  elastic 
axis  and  the  engine.  A calculated  time  history  of  the  acceler- 
ation at  the  accelerometer  location  is  given  in  figure  14;  this 
time  history  is  based  on  linear  interpolation  between  the 
computed  accelerations  at  the  engine  and  the  clastic  axis  at 
this  station  in  accordance  with  the  proportionate  distances 
involved.  The  figure  also  includes  the  recorded  maximum 
acceleration  at  this  point  and  shows  good  agreement  of  the 
computed  acceleration  therewith. 
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Acceleration  time  histories  for  the  engines,  computed  on 
tile  basis  of  a rigid  structure,  are  included  in  figure  14.  The 
maximum  accelerations  computed  on  this  basis  do  not  agree 
with  the  structural  failures  which  occurred.  Furthermore, 
the  fact  that  the  curves  computed  on  the  basis  of  a rigid 
structure  reach  a maximum  at  the  same  instant  of  time  is  in 
strong  disagreement  with  experiment,  which  in  this  respect 
is  in  approximate  agreement  with  the  computations  for  ah 
elastic  structure. 

Difference  between  the  ratio  f*/f(  for  the  inboard-engine 
and  outboard-engine  stations  is  primarily  responsible  for 
difference  in  the  calculated  response  of  the  engines.  Most 
of  the  difference  in  this  ratio  for  the  two  stations  is  not  due 
to  difference  in  the  sprung-engine  frequency  but  is  due  to 
the  greatly  different  time  to  reach  maximum  acceleration  ft 
of  the  elastic  axis.  Agreement  with  experiment  of  the 
response  calculated  at  these  stations  provides  indirect  con- 
firmation of  the  acceleration  time  histories  predicted  for  the 
elastic  axis  at  these  locations  by  the  normal-mode  method. 

For  the  impact,  experimental  data  are  not  available  for 
checking  the  tip  acceleration,  which  is  predicted  on  the 
basis  of  the  normal-mode  method,  but  the  initial  downward 
acceleration  and  the  12 g maximum  acceleration  shown  in 


figure  13  for  this  station  are  in  general  agreement  with 
results  recorded  in  severe  impacts  of  other  airplanes.  Agree- 
ment of  the  computed  hull  acceleration  with  experiment  has 
already  been  shown;  in  an  indirect  manner  all  the  curves  in 
figure  13  exhibit  satisfactory  agreement  with  available  ex- 
perimental data.  Although  the  response  of  the  engines  is 
different  from  the  response  assumed  in  calculating  these 
curves,  it  appears  that  in  practical  use  the  two-mass  solu- 
tions given  herein  can  be  interpreted  on  the  basis  of  the 
normal-mode  method  to  obtain  both  the  response  of  the 
' elastic  axis  and  the  hydrodynamic  force.  Further,  it  appears 
that  in  practical  problems  modification  of  these  two  results 
is  not  required  when  accelerations,  of  the  engine  different 
from  the  accelerations  predicted  by  the  normal-mode  method 
are  determined;  the  merit  of  this  statement'  should  be 
independent  of  whether  such  modification  is  made  by  the 
method  used  herein  or  by  another  method. 

In  the  foregoing  comparison  agreement  of  calculation  with 
experiment  is  obtained  without  consideration  of  the  response 
of  modes  higher  than  the  fundamental  mode.  If,  when 
more  data  are  available,  it  is  shown  that  the  response  of  the 
higher  modes  can  be  determined  by  treatment  parallel  to 
that  given  the  fundamental  mode,  the  two-mass  solution 
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given  herein  can  be  used  to  predict  their  response  by  selecting 
a solution,  for  a mass  ratio  and  natural  frequency  reprcsenta- 
tive  of  the  higher  mode,  which  has  a nodal  acceleration  curve 
of  approximately  the  same  shape  as  the  hydrodynamic-force 
curve. 

CONCLUSIONS 

Theoretical  solution  of  hydrodynamic  impact  of  a hull 
mass  connected  by  a spring  to  an  upper  mass,  the  results  of 
calculations  for  wide  ranges  of  mass  ratio  and  natural  fre- 
quency, and  the  use  of  these  results  in  a comparison  of 
theoretical  data  with  test  data  for  a flight-tost  landing 
impact  indicated  that:  * 

1.  In  flying-boat  impact  the  effect  of  the  structural 
response  on  the  hydrodynamic  force  might  be  substantial, 
the  shape  of  the  forco  time  history  might  bo  considerabl}' 
changed,  and  the  maximum  hydrodynamic  force  might  be 
either  reduced  or  increased. 

2.  The  greatest  reduction  in  hydrodynamic  force  occurred 
for  the  condition  of  large  mass  ratio  and  1owt  value  of  spring 
constant. 


3.  The  normal-mode  method  was  a practical  means  for 
determining  the  equivalent  two-mass  system  which  repre- 
sented the  major  elastic  action  of  the  airframe,  for  predicting 
the  effect  of  this  action  on  the  hydrodynamic  force,  and  for 
approximating  accelerations  along  the  elastic  axis  of  the 
wing. 

4.  The  acceleration  of  engines  contained  in  nacelles  for- 
ward of  the  wing  could  not  he  computed  on  the  basis  of 
coupling  between  torsion  and  bending  as  in  the  fundamental 
mode,  but  simple  treatment  of  their  response  to"  the  com- 
bined translation  and  oscillation  of  the  elastic  axis  gave 
agreements  with  experimental  accelerations  and  gave  an 
explanation  of  partial  failure  of  the  inboard-engine  mounts 
during  impact. 


Langley  Memorial  Aeronautical  Laboratory, 
National  Advisory  Committee  for  Af.ronautics,  . 
Langley  Field,  Va.,  March  17,  1947. 
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APPENDIX  A 

MATHEMATICAL  EQUATIONS  AND  METHOD  OF  SOLUTION 
EQUATIONS  OF  MOTION  ( COMPUTING  DIRECTIONS 


The  following  equations  of  motion,  which  were  derived 
from  equation  (30)  in  reference  2,  are  for  fixed-trim  impacts 
of  a rigid  prismatic  float  connected  by  a massless  spring  to 
a rigid  upper  mass: 

Acceleration  oij  lower  mass  normal  to  water  surface  in  feet 
per  second  per  second 

-y^\zAyIJ\yL+K1cosT)t+^-s  [^+ 

(ms+mI,)(yL—y^t)+2K1Acosr  £ yL3dt+ 

(AI) 

Acceleration  of  sprung  mass  normal  to  water  surface  in 
feet  per  second  per  second 

Vs=~r  l(AyL3 + mL)  y L -f 3 A yL-  (yL-\-Kt  cos  r)*]  (A2) 

771  g 

Acceleration  of  nodal  point  normal  to  water  surface  in  feet 
per  second  per  second 


- __rnLyL+msys 
mL+m3 


(A3) 


The  spring  constant  in  pounds  per  foot  of  deflection  can  lie 
expressed  by  the  equation  1 


j.r_±T-(mLma)jAi  ■ 

mL+ms 


(A4) 


where 

A’1=ii0sin  r — yLi 


0 cos  T 


A = 0.82 


(35-OX 


YV  i 

tan  - \ 

/ rp  \ 

)\l 

2 tan  j3  / 

\6  sinr  cos*r/ 

yL  velocity  of  lower  mass  normal  to  water  surface,  ft/sec 
yL  draft  normal  to  water  surface,  ft 
/3  angle  of  dead  rise,  radians 

r angle  of  trim,  deg 

yLl)  initial  velocity  normal  to  water  surface,  ft/sec 
XcQ  initial  velocity  parallel  to  water  surface)  ft/sec 
rn3  sprung  mass,  slugs 
mL  lower  mass,  slugs 
/,  natural  bending  frequency,  cycles/sec 
p mass  density  of  water,  slugs/cu  ft 


A sample  data  sheet  is  given  as  table  I.  In  this  table  the 
numbers  in  circles  refer  to  rows;  the  circled  numbers  under 
the  row  headings  refer  to  computed  values  to  be  used  for 
the  computations.  In  the  first  column  the  time  is  equal  to' 
zero,  in  the  second  column  the  time  is  At,  and  in  successive 
columns  the  time  is  2 At,  3 Af,  and  so  forth.  Each  row  is 
computed  in  sequence  for  any  given  column  before  any  row 
is  computed  for  the  next  column,  except  for  the  first  column, 

TABLE  I.— DATA  SHEET— GENERAL  TERMS 


I Row- 


Row  heading 


Column 

1 t 

r-  -r: 

1:2,3 

A 

l 1 

' "111 

! L 


1 @ 


@ 


Tune,  sec 

assumed 

t@,XA0-K*>»=IfL,ft 

t=Vi} 

=Vl* 

<+Kl  cos  r)J@X3n 
®XT 

(nti.+  ms)  (®— ifi^X®) 

2JVI.4  COS  rX@ 

’~(®+®,)fj+@P=J’tin1A  * 

3.4AV  cos*  rX@  . 

®-f®+@+@ 

(f.x®5+® 

(AX®)  + ffli 
— ft/sec5 

(@+@p)f 

@ + @p 

vi^A® 

®+f 

vl  assumed  for  next  interval 

r_-L(©-@)i 

L m&  J . 

32J units 

325=#!.,  ? units 
("itXg)-f(maX8). 

771 lA  77 Is 


S units 
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the  tune  for  a given  mass  to  reach  a maximum  acceleration.) 

The  number  of  significant  figures  to  be  used  should  ho 
chosen  on  the  basis  of  the  computing  equipment  available 
and  the  accuracy  desired.  Comparison  of  a solution  com- 
puted with  four  significant  figures  with  results  which  had 
been  obtained  with  six  significant  figures  gave  a difference  of 
about  1 percent  in  the  maximum  acceleration. 

As  the  computations  proceed,  the  lower-mass  acceleration 
(row  ©),  the  sprung-mass  acceleration  (row  @),  and  tho 
nodal-point  acceleration  (row  ©),  all  in  g units,  should  be 
plotted  against  time,  (row  ©)  in  seconds.*  Each  time  a 
column  is  completed  the  new  points  should  be  added  to  the 
plot.  This  plot-  is  the  only  brief  method  of  checking  on  the 
accuracy  of  the  computation  until  sufficient  solutions  are 
obtained  to  permit  cross-plotting.  If  the  points  do  not  lie 
on  a smooth  curve,  an  error  lias  been  made  in  tho  computa- 
tions. If  an  error  is  made,  every  value  computed  thereafter 
contains  the  error  and  therefore  great  care  must  be  taken. 

The  number  of  columns  required  for  the  computations 
depends  on  tire  number  of  columns  required  for  rows  © and 
© to  pass  their  respective  maximums.  At  least  three 
columns  should  be  computed  beyond  the  column  in  which 
row  @ reaches  its  maximum  value. 

TABLE  II.— SAMPLE  COMPUTATION 


[Comtanta  are  from  appendix  A] 


— . 

Row. 

• 

Row  heading 

— ••••  aJg  ;■  a.- 

Column 

1 

2 

3 

4 

5 

6 

. l 

7 

8 

Qn 

@ 

@ 

Time  At,  sec  (0.005)" 

®v — — - 

(@PX0.005)+®P 

(®+ 4.23736“)  X 401 .757 

®X  33. 4798 

1240.991:®—  (20.6673  X ©)  ] 

F(®+®*)0.0025]+@J, 

1 134.92  X® - 

f(®+fflI,)0.0025]+@, 

K®+@p)0.0025H-@b 

7213.62X® - 

®+©+@+® 

(150.533X@)  + ® 

(133.9I9X®) +525.776 



(@+@^0.0025 

20^6^+® i 

®+f-  ....... ■ ■ 

@+f ----- — 

r-0.  00l398(®—@)  1/32.2. 

@/32.2 

(625.776  X@)  + (715.217X®) 
1240.99 

.0 

20.  6673 
0 

...  - 0 

0 

. 0 
0 
0 
0 
0 
0 

- 0 
0 

-0 

, ' .0 
0 

- 0 

,.D 

0 

20.  6673 
20.  6073 

20.  6673 

J) 

-.0 

- • • -0 

6.015 
20.  5287 
.0.  309380 
0.  095716 
0.  029613 
6.  003162 
23586.5 
0.  306742 
-0.  781824 
0.  000124 
0.  140730 
0.  000506 
0.  000002 
0.  014427 
-0.  319925 
2353$.  $ 
529.  7i2 
-44.  4335 
-0.  161415 
— 0..  237056 
20.  4302 

20.  3455 

20.  268£ 

-0.  002tfij3 
-1.  37992 

-0.  585844 

0.  020 
20.  2688 
0.  411128 
0.  169026 
0.  069491 
0.  028570 
40782.0 
0.  956518 
-2.  75252 
0.  000373 
0.  422190 
0.  001168 
0.  000006 
0.  043282 
-0.  133053 
40581.  7 
535.  082 
-75.  8420 
-0.  300690 
-0,  537746 
20.  1296 

19.  9793 

’ 19.  8290 

-0.  008696 
-2.  35534 

-1.  00291 

il»V  1 • — .-r  - V_J3 

0.  025 
19.  8290 
0.  511025 
0.  2611147 
0.  133453 
0.  0G8198 
60767.  4 
2.  28326 
-7.  02152 
0.  000879 
0.  997595 
0.  002243 
0.  000015 
0.  108204 
-3.  63246 
60220.  6 
543.  648 
-110.  771 

-0.  466533 
-1.  00428 
19.  6630 

19.  4297 

19.  1964 

-0.023740 
-8  44009 

-1.  47110 

ia 

0.  030 
19.  1964 
0.  608174 
■ 0.  369876 
0.  224949 
0.  136808  ' 
81602.  9 
4.  58030 
-14.  8996 
a 001775 
2.  01448 
0.  003821, 
0.  000030 
0.  216409 
-7.  88831 
80415.  5 
555.  901 
-144.  658 
-0.  638573 
-1.  64285 
19.  0245 

18.  7052 

18  3859 

-0.051552 
-4.  49248 

-1.  93306, 

0.  035 
18  3859 
0.  701700 
0.  492383 
0.  345505 
0.  242441 
101246.  0 
8 11688 
—26.  8749 
0.  003201 
3.  63288 
0.  005977 
0.  000054 
0.  389525 
-14.  7356 
99028  0 
572.  040 
-173.  112 
-ft  794425 
-2.  43728 
18  2300 

17.  8328 

. 17.  4355 

-0.  090297 
-6.  37616 

-2.  33323 

In  the  first  column  the  assumed  value  of  yL  is  the  given  value 
of  ifi^.  All  other  values  are  zero  except  the  values  for  tows 
@,  and  @ which  in  this  case  are  also  y^.  Lower-case  p 
as  a subscript  on  a row  number  refers  to  the  computed  value 
in  the  indicated  row  of  the  preceding  column. 

Each  row  heading  indicates  tho  general  operation  to  be 
performed.  When  the  data  sheet  is  set  up,  numerical  values 
should  be  substituted  in  the  row  headings  for  the  expressions 
that  are  constant  for  a specific  case.  The  constants  for  the 
sample  computation,  given  in  table  II,  are  as  follows: 


K - 107,  664 

Ki. 4.24317 

.4 -133.  919. 

ft,  deg — -=-  > 22.-S 

t,  deg — — . ..  .<  -3 

ft/uec 20.  6673 

ft/eec -> — -----  82.  1685 

ms ..  715.217 

mL — — . 525.  776 

/«,  cyeles/sec ■ 3 

p,  slugs/cu  ft ... - , L 938 


For  most  solutions  a time  increment  At  of  0.005  second  is 
satisfactory.  (The  value  of  At  may  he  varied  by  considering 
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APPENDIX  B 

RELATIONS  BETWEEN  TWO-MASS  SYSTEM  AND  REPRESENTED  STRUCTURAL  MODE 


Tlie  sum  of  the  masses  mL  and  ms  shown  In  figure  1 must 
equal  the  gross  mass  of  the  represented  airplane  in  order  to 
obtain  the  proper  nodal  acceleration.  For  the  hypothetical 
limit  condition  in  which  the  wing  mass  is  concentrated  at  a 
single  point  in  each  semispan,  mL  is  the  actual  hull  mass  and 
ms  is  the  actual  wing  mass.  In  order  to  take  into  account 
the  more  complex  nature  of  the  structural  action  for  a par- 
ticular mode,  the  determination  of  the  ratio  of  these  masses 
is  necessary  so  that  the  vibrational  energy  of  the  simplified 
and  represented  systems  are  equal  for  the  same  vibrational 
amplitudes  of  mL  and  the  actual  hull  or  point  of  force 
application. 

On  the  basis  of  the  theory  of  vibrations,  for  which  equa- 
tions are  included  in  reference  3,  the  vibration  energy  E of 
the  two-mass  system  is  given  by  the  equation 

^5=2  2 ms<Psi  (Bl) 

where 


Equality  of  the  vibrational  amplitude  of  the  lower  mass  of 
the  simplified  system  to  the  vibrational  amplitude  of  the 
hull  or  fuselage  of  the  flying  boat  or  airplane  relative  to  the 
nodal  point  of  the  represented  mode  is  expressed  by  the 
equation 

<PL=<Pk  (B5) 

where 

<pi  deflection  of  hull  or  fuselage  of  flying  boat  or  airplane 
relative  to  nodal  point  of  represented  mode 

The  requirement  of  equal  energy  of  the  simplified  and 
represented  systems  for  the  condition  expressed  by’  equation 
(B5)  gives  combination  of  equatipns  (B3),  (B4),  and  (B5) 
to  obtain  the  following  equation  for  the  mass  ratio  of  the 
two-mass  system: 

®s= , flgg. 


•fL  vibrational  amplitude  of  mL  relative  to  nodal  point  of 
system 

vibrational  amplitude  of  ms  relative  to  nodal  point  of 
system 

to  natural  frequency 

Since  for  the  two-mass  system  the  node  is  at  the  center  of 
gravity, 

= Wl,g^s  (32) 


If  equation  (Bl)  and  equation  (B2)  are  combined,  and  since 
the  total  mass  m is  equal  to  the  sum  of  the  masses  ms  and 
mL,  the  following  equation  can  be  obtained  for  the  vibrational 
energy  E of  the  two-mass  system: 


E 1 mmh  5 
«*  2 ms  <?L 


(B3) 


The  vibrational  energy  of  the  represented  mode  is  a function 
of  the  spanwise  mass  distribution  and  mode  shape.  On  the 
basis  that  the  semispan  of  the  airplane  is  divided  into  j sec- 
tions or  stations,  the  vibrational  energy  of  the  mode  can  be 
written  as  follows: 


Computation  of  the  natural  frequency  of  wing  modes  has 
received  a great  deal  of  attention  in  connection  with  study 
of  wing  Sutter  and  need  not  be  treated  herein.  Incidental 
to  calculation  of  the  natural  frequency,  a mode  shape  is 
attained  which,  together  with  knowledge  of  the  mass,  distri- 
bution, permits  use  of  equation  (B6).  In  cases  in  which  the 
wing  has  been  constructed,  the  mode  shape  and  natural 
frequency  may  be  determined  experimentally.  Equation 
(A4)  in  appendix  A of  the  present  report  permits  computa- 
tion of  the  spring  constant  which  for  a given  mass  ratio  of 
the  simplified  system  gives  the  required  natural  frequency. 

After  the  accelerations  of  the  two  masses  of  the  simplified 
system  have  been  computed,  equation  (A3)  of  appendix  A 
fixes  the  magnitude  of  the  nodal  acceleration.  The  difference 
between  the  nodal  acceleration  and  the  hull,  fuselage,  or 
float  acceleration  can  he  taken  as  a measure  of  the  oscillatory’ 
acceleration.  On  the  basis  that  the  structure  deflects  in  the 
mode  used  in  determining  the  equivalent  two-mass  system, 
the  acceleration  at  any  point  is  given  by  the  equation 

Ttp=nt+n,—  - - (B7) 

<Pk 


^=4  (m1<p11+m2<fhi+m3a>it+  . . . + (B4) 

03“  Z 

where 

m}  mass  of  j’tli  spanwise  section;  value  for  semispan  doubled 
to  represent  the  entire  span 

tp}  deflection  of  mass  at  j’tk  spanw'ise  section  relative  to 
nodal  point  (h}+xaf) 

■ hj  deflection  of  elastic  axis  at  jth  spanwise  deflection  rela- 
tive to  nodal  point 

x chordwise  distance  from  elastic  axis  to  effective  mas3 
center 

a.)  torsional  deflection  at  j’th  spanwise  station 


The  foregoing  equations,  with  consideration  of  rotatory 
inertia  and  energy,  may  be  applied  to  the  case  of  a tip  float 
attached  to  a flexible  wing  if  the  stiffness  and  mass  distribu- 
tion are  known  and  if  a manner  of  structural  deflection  is 
assumed. 
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